计算物理 ›› 2021, Vol. 38 ›› Issue (1): 5-15.DOI: 10. 19596/ j. cnki. 1001-246x. 8237

• 研究论文 • 上一篇    下一篇

磁瑞利-泰勒不稳定性非线性演化数值模拟

曹启伟1,2, 肖德龙1, 杨显俊1, 王建国1   

  1. 1. 北京应用物理与计算数学研究所, 北京 100094;
    2.
    中国工程物理研究院研究生院, 北京 100083

  • 收稿日期:2020-05-29 修回日期:2020-06-07 出版日期:2021-01-25 发布日期:2021-01-25
  • 通讯作者: 肖德龙,E-mail:xiao_delong@iapcm.ac.cn
  • 作者简介:曹启伟(1995-),男,硕士生,从事磁惯性约束聚变研究,E-mail:caoqw17@163.com
  • 基金资助:

    国家自然科学基金(11775032)、国家重点研发计划(2017YFE0301801)资助项目

Numerical Simulation on Nonlinear Evolution of Magneto-Rayleigh-Taylor Instability

CAO Qiwei1,2, XIAO Delong1, YANG Xianjun1, WANG Jianguo1   

  1. 1. Institute of Applied Physics and Computational Mathematics, Beijing 100094, China;
    2. Graduate School of China Academy of Engineering Physics, Beijing 100083, China


  • Received:2020-05-29 Revised:2020-06-07 Online:2021-01-25 Published:2021-01-25
  • Contact: Xiao Delong,E-mail:xiao_delong@iapcm.ac.cn

摘要:

    基于理想磁流体力学模型,采用5WENO有限差分格式结合磁场伪散度清除投影算法,发展二维数值模拟程序,利用线性理论对程序进行检验,在此基础上研究磁瑞利-泰勒(MRT)不稳定性非线性演化物理过程和发展规律。结果表明:垂直于加速度方向磁场在线性和非线性阶段对MRT扰动发展都有较强的制稳作用,而平行于加速度方向磁场在线性阶段对MRT扰动发展影响较小,但在MRT扰动非线性发展阶段会对由于流场变化产生的K-H不稳定性有显著的制稳作用,从而抑制总的扰动发展。在MRT非线性发展过程中,由于波模耦合相互作用,单模演化会产生高次谐波;双模演化会产生除基模和谐波以外的其他波模,其扰动幅度甚至超过基模;多模演化过程中发生逆级联过程,即短波长扰动逐渐向长波长扰动发展。


关键词:

磁流体力学, WENO有限差分格式, 磁瑞利-泰勒不稳定性, 线性稳定性, 波模耦合

Abstract:

A 2D numerical simulation program is developed for nonlinear evolution of magneto-Rayleigh-Taylor instability(MRTI). Based on an ideal magnetohydrodynamic model, a 5th WENO scheme is employed to calculate MHD equations and a projection scheme is used to clean spurious magnetic field divergence. The method is validated with MRTI compared with theory of linear stability. It is shown that the vertical magnetic field suppresses evidently MRTI in both linear and nonlinear stages. The parallel magnetic field has weak effect on MRTI in linear stage. However, it reduces remarkably Kelvin-Helmholtz instability and overall instability in nonlinear regime. High order harmonic occurs in the nonlinear evolution of single mode, and many modes except fundamental modes and harmonic occur in the evolution of double modes. In MRT evolution of multi-mode seeds inverse cascade occur, which means that perturbations evolve from short wavelength to long wavelength in the nonlinear stage.

Key words: Magnetohydrodynamics, WENO finite difference scheme, magneto-Rayleigh-Taylor instability, linear stability, mode coupling